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Abstract 

We consider the one-dimensional anisotropic XY model in the continuum limit. Stability analysis 
of its Bloch wall solution is hindered by the nondiagonality of the associated linearised operator 
and the hessian of energy. We circumvent this difficulty by showing that the energy admits a 
Bogomolnyi bound in elliptic coordinates and that the Bloch wall saturates it — that is, the Bloch 
wall renders the energy minimum. Our analysis provides a simple but nontrivial application of the 
BPS (Bogomolnyi - Prasad - Sommerfield) construction in one dimension, where its use is often 
believed to be limited to reproducing results obtainable by other means. 
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I. INTRODUCTION 



Although the Bogomolnyi decomposition is an indispensable tool for the study of higher- 
dimensional field theories it has seldom been used in one dimension. It does allow 
one to reduce the order and find topological solitons of one-component models such as 
the sine-Gordon and 4 -theory, but these can be obtained with less effort simply by using 
an integrating factor. It is of more value for multicomponent and lattice jsj] systems, 
but all models that benefited from its use had to be specifically designed to admit such a 
decomposition. 

In this paper, we apply the Bogomolnyi construction to a system that has been studied 
for more than forty years: the anisotropic XY model. It has long been believed that the 
Bloch-wall solutions of this model are stable; this fact has been demonstrated numerically 
but never proven analytically. Here, by transforming to elliptic coordinates, we find the 
Bogomolnyi bound for the energy and show that the Bloch wall minimizes the energy in 
the corresponding topological sector; this proves its stability. (Note that the Bogomolnyi 
construction cannot be carried out in the original coordinates.) 

An outline of the paper is as follows. In the next section we introduce the (station- 
ary) anisotropic XY model and its domain-wall solutions. We also describe three different 
time-dependent extensions of the model which correspond to three different types of be- 
haviour, viz., dissipative, conservative relativistic and conservative nonrelativistic dynamics. 
In section HTT1 we show that because of the nondiagonality of the linearised operator and the 
second variation of energy, the standard methods are inadequate for proving the stability of 
the Bloch wall — in each of the three cases. Section [IV] contains the main result of this pa- 
per: the proof of the energy minimization by the Bloch wall. Finally, section IVl summarises 
the results of our work. 



II. THE MODEL 



A. Anisotropic easy-axis ferromagnet near the Curie temperature 



The anisotropic XY model was originally introduced to describe domain walls in an easy- 
axis ferromagnet near the Curie point 4, s|. In the study of domain walls, the magnetization 



is assumed to vary only in x direction and the magnetization vector M is taken to lie in the 



yz plane. (With this notation, the system should have been called the Y Z model but we 
keep the traditional name to avoid confusion.) In the continuum limit, the model is defined 

fin 

by its free energy expansion |4|, |5j : 



E 



l -(d x M) 2 - (1 + h)M 2 + ^M 4 + 2/iM 2 + £ 



dx. 



(1) 



The first term in ([I]) is the exchange energy, which is minimized when M = const. The 
combination of the second and third terms makes the |M| ^ ground state energetically 
preferable to the state with |M| = 0. Anisotropy is caused by the term 2/iM 2 ; the parameter 
h is assumed positive, making the z direction the easy-axis. Finally, the coefficient (1 + h) 
is introduced to simplify the subsequent formulas and a constant £ has been added to the 
integrand to ensure that the free energy is finite. 

Defining x/j = M y + iM z , we can recast the free energy in the form 



E 



-M 2 + 



\ 2 + ^ 2 + r 2 )+s 



dx. 



(2) 



This is the form that we shall be working with in this paper. 



B. Anisotropic easy-plane ferromagnet in external field 

The same free energy expansion (T5]) describes one more, unrelated, magnetic system: a 
weakly anisotropic easy-plane ferromagnet in an external field perpendicular to the easy 



plane 



6|. In terms of the magnetization vector M, the free energy of this system is given by 



E 



\{d x M) 2 + \Ml + ^Ml - HM Z + £ 



dx, (3) 



where H > is an applied field, the anisotropy parameters (3 and e/3 are positive, and 
e < 1. In this model, the xy plane is the easy plane, with a weak anisotropy favouring the 
y direction. Far from the Curie point, the magnetization is hard, i.e. the magnitude of the 
magnetisation vector is constant: M 2 = Mq = const. [For this reason equation Q has no 
terms in M 2 and M 4 .] 

We assume that the the magnitude of the external magnetic field is close to f3M$: H = 
(3 Mo — eq, where q is a quantity of order 1. The (degenerate) ground states will then have 
the magnetization vector almost parallel to the field: 



The deviation from the reference vector M = (0, 0, Mo) pointing in the direction of the 
field can be characterised using the complex variable ip = (f3/2es) 1 ^ 2 (M x — iM y ), where 
s = qM - 15 Mil 2. The M z component of the magnetisation is expressible using M = 
M5 2 = const: 



(3M> 

Transforming to x = (es/2aMl) 1 / 2 x and keeping terms up to e 2 , the free energy (^j reduces 
to equation (j2j) where h = /3Mq/2s and we have dropped the tilde over x. Unlike the 
anisotropic XY model considered in the previous subsection, the domain walls in the present 
system interpolate between the nearly-parallel ground states M^ + ^ and M>'. It is fitting 
to note here that due to the relative weakness of magnetic anisotropy, the field H pa /3M 
will generally be well within the range of modern experiments. For example, the martensite 
phase of NiMnGa is a weakly anisotropic easy-plane ferromagnet with (3 ~ 2 x 10 _6 Tm/A 
and M Q pa 5 x 10 5 A/m |7j. Therefore, when subjected to a magnetic field H pa IT, it will 
be described by equation (j2J). 

Finally, we note that the expansion ([I])-© arises in yet another magnetic context, viz. a 
strongly anisotropic ferromagnet with spin S = 1 



C. Bloch and Ising walls 

The free energy is extremised by solutions of the following stationary equation: 

\ij>sx-\ij>\ 2 il> + ip-ty* = 0. (4) 

There are two soliton, or kink, solutions available in literature; each describes an interface 
between two ferromagnetic domains. One is commonly known as the Ising, or Neel, wall 



2, q, y, 



10|: 



ipi(x) = iAtanh(Ax), (5) 



where A = yl + h. The second kink solution has the form 

^Pb(x) = %A tanh(Bx) ± Csech(Bx), (6) 

where B = \^4h and C = \/T — 3h , 121. This solution is referred to as the Bloch wall. 
The two signs in front of the real part in ([6]) distinguish Bloch walls of opposite chiralities. 
The Ising wall exists for all positive h whereas the Bloch wall exists only for h < g. 



The Bloch and Ising wall have the same asymptotic behaviour: ipB,i{x) — > ±iA asi-> 
±00. This determines the constant So that is added to the integrand in (j2J) to ensure the 
finiteness of the free energy: £ = A A /2. The Bloch and Ising walls have energies 

E B = 2B(A 2 - B 2 /3) (7) 

and Ej = |^4 3 , respectively. It is easily verified that the energy of the Ising wall is greater 
than the energy of the Bloch wall for all h < |. 

In addition to the Bloch and Ising walls, equation (jlj) has a family of nontopological 
solitons, available in explicit form and describing Bloch-Ising bound states [13]. These will 
not be considered here. 



D. Three types of dynamical behaviour 



The dynamics of the domain walls is governed by one of three possible time-dependent 
extensions of equation (jlj). In the case when the free energy (j2J) is used to model the 
ferromagnet near the Curie point, the evolution of the field ip is dissipative and governed by 



the Ginsburg-Landau equation 



ip t = ip 



hip*. 



(8) 



On the other hand, in the case of the ferromagnet with spin 5 = 1, the field ip satisfies a 
relativistically-invariant equation 



^-^-2^ + 21^ + 2^* 



0. 



(9) 



Originally, this equation was introduced by Montonen 11[ in a different context — as an 
exactly solvable special case of Rajaraman and Weinberg's bag model 9J. Independently, 
Sarker, Trullinger and Bishop [12| proposed it as an interesting interpolate between the 
sine-Gordon and the 4 theories. Accordingly, the Klein-Gordon equation is commonly 
known as the Montonen-Sarker-Trullinger-Bishop (MSTB) model. 

Finally, the magnetisation vector of an anisotropic easy-plane ferromagnet in an external 
field satisfies the Landau-Lifshitz equation. The perturbation procedure described in section 
III Bl reduces it j(3] to the parametrically driven nonlinear Schrodinger equation, 



iipt + - v#r + v> = w- 



(10) 



The stability of the Bloch or Ising wall depends on which of the three equations ©, © 
or fflOl) governs the evolution of ip. Consequently, the three cases — the Ginsburg-Landau, 
the Klein-Gordon and the nonlinear Schrodinger — need to be considered separately. 



III. APPROACHES TO STABILITY 

In this section we describe the standard methods of stability analysis and explain why 
they all fail in the case of the Bloch wall — for each of the three types of evolution. 

A. The Ginsburg-Landau and relativistic dynamics 

We start with the Ginsburg-Landau equation, equation ([8]), and linearise it about the 
stationary solution ipo{x), which can be either the Bloch or Ising wall. Decomposing the 
small perturbation 6ip(x, t) into its real and imaginary parts, 5ip(x, t) = 81Z(x, t) + i52(x, t), 
and letting 

6H(x, t) = 5TZ(x)e x \ 5I(x, t) = 51(x)e x \ 
yields an eigenvalue problem 





-n \=x , (ii) 



where 7i is a self-adjoint operator 




w= l Idl-i + h + ml + i^ 2^ Xo | (12) 



l d 2_ 1 _ h + 1Z 2 + 3I 2 



In (112p . TZq and 1q stand for the real and imaginary part of the stationary solution ipo, 
respectively. The solution ip will be unstable if the operator — 7i has at least one positive 
eigenvalue A, and stable otherwise. 

In the case of the MSTB model [equation (jUJ)], the linearisation about ipo produces the 
same eigenvalue problem (TTTI) . with the same operator ( Tl2l . where one just needs to replace 
A with A 2 /2. Here we have the same stability criterion as in the Ginsburg-Landau case: the 
solution ip will be unstable if — TC has at least one positive eigenvalue A 2 /2. 

An alternative approach (which still leads to the same criterion, though) is based on 
considering the Lyapunov functional (see e.g. [15|). The Ginsburg-Landau equation (jSJ) can 



be written as ipt = —SE/dip*, where E is the functional ([2]). Hence this functional satisfies 

E t = -2 / \ij t \ 2 dx, 



and so we have E t < unless ip is a static solution, ip t — 0. Now if we could prove that 
£/[^] > -E^o] f° r an V* m some neighbourhood of ?/>o, the functional would be the Lyapunov 
functional for this solution and hence ipo would be proven to be stable. 

Here we need to make a standard remark on the translational invariance of equations (jSJ), 
(Q and (jTUl) . The domain wall centred at the origin has the same energy as the wall centred 
at any other point xq and therefore can never be an isolated minimum of E. However, 
a mere translation of the wall from the origin to the point xq does not imply instability. 
Therefore, it is physically reasonable to group all configurations obtained from a given ip{x) 
by translations if)(x) — > ij){x — x ) with — oo < x < oo, into equivalence classes. Unlike 
the wall with any particular xq, the equivalence class consisting of all translated walls can 
be an isolated minimum of the energy — defined on the corresponding quotient manifold. 
If the perturbations of the wall are assumed to be infinitesimal, the quotient is a linear 
subspace; namely, it is the quotient of the space of all infinitesimal perturbations of the wall 
by the subspace spanned by its translation mode (d x TZo,d x X ). This quotient space can be 
conveniently characterised by the orthogonality constraint 



(d x K ,dM\ \dx = 0. (13) 




In what follows, the subspace of perturbations defined by the constraint ([TBI will be denoted 
S. 

To check whether ipo renders E[ip] a minimum in <S, the functional is expanded about the 
stationary point ip Q \ 

E[ij] =E + ~5 2 E + .... 
Here ip = ipo + 5TZ + i5I, and the second variation has the form 

I«5 2 £ = J (6K, 51) H ( 6 ™ j dx, (14) 

where the hessian Ti coincides with the linearised operator (112p . The solution ipo will min- 
imise E[ip] in S provided TC has no negative or zero eigenvalues other than the one associated 
with the translation mode. 



7 



A Lyapunov functional can also be used in the case of the relativistic dynamics © (see 
e.g. 



161]) . Here, as a candidate functional one considers the total energy 

£totai[^t] = \f \A\ 2 dx + E[ij>], (15) 

where E[i/j] is as in §2$. The energy is conserved, and hence E tot3i \ will define the Lyapunov 
functional for the solution ipo if it renders this functional a minimum in S. Since the first 
term in (fT5l) is minimised by any static configuration, it is sufficient to check whether ipo{x) 
renders the functional E[ip] a minimum. Consequently, the minimisation problem reduces 
to the eigenvalue problem (TTTT) . 

Thus, we have the same stability criterion in the case of the Ginsburg-Landau and rela- 
tivistic dynamics — one has to prove that the solution minimises the functional ([2]) under 
the constraint ( |T3l) . or, equivalently, show that the operator 7i has no negative or zero 
eigenvalues other than the translation mode. 

In the case of the Ising wall ([5]), we have 1Zq(x) = and the operator (|T2|) is diagonal. Its 
eigenvalues can be readily found and the above criterion easily implemented. Thus, it was 



17 



18| that the Ising wall is stable for h > | and unstable for h < ~. (See 



shown in 

{19I ] for the generalisation to a nonvariational case.) 

In the case of the Bloch wall, on the other hand, the operator 7i is nondiagonal. This 
makes it impossible to determine the sign of the lowest eigenvalue of 7i using standard 
analytical methods. Consequently, previous studies had to resort to semi-intuitive and 



numerical arguments. In particular, it was noted that the energy of the Bloch wal 
than that of the Ising wall and suggested that the Bloch wall should be stable [ll 



is 



12 



owcr 



20]. 



This conjecture was supported by results of direct numerical simulations of equation ([9]) 



2l| and perturbation theory for small C (C — \/\ — 3h) ]8]. However no analytical proof, 



applicable for all parameter values, has been given so far. 



B. The Schrodinger dynamics 

Finally, we discuss the parametrically driven nonlinear Schrodinger equation, equation 
TUT) . In this case, the linearisation about x/jq produces a symplectic eigenvalue problem 

J=|°~^|. (16) 





where Tt is as in ([12]) . The product operator J~ l TL is non-self- adjoint and hence its eigen- 
values can be complex. The solution ipo will be unstable if the operator J~ l 7i has at least 
one eigenvalue A with positive real part. 

The nonlinear Schrodinger equation ffTUj) conserves energy which is given by the integral 
(J2J); hence the energy is a potential Lyapunov functional for equation (1101) . One can therefore 
try to establish the stability of ipo by proving that it renders the energy minimum under the 
constraint ([131) ; this happens when the operator 7i in the second variation fTHl) does not have 
negative or zero eigenvalues other than the translation mode. Note that the criteria based on 
the linearisation and the energy minimality appeal to eigenvalues of different operators here, 
J~ X H and 7i, respectively. However, it is not difficult to show that the positive definiteness 
of 7i implies that J^ l 7i does not have eigenvalues with positive real part. We include a 
proof of this simple fact in the Appendix. 

As we have pointed out in the previous subsection, the operator 7i associated with the 
Ising wall is diagonal. Making use of this property, it was proved in [6| that J~ 1 H. does not 
have eigenvalues with positive real part and the Ising wall is stable for all h > 0. On the 
other hand, in the case of the Bloch wall, the operator 7i in ( TTTf) and ( fl6l) is nondiagonal. 
This prevents the determination of the sign of the lowest eigenvalue of 7i, or testing the 
existence of unstable eigenvalues of J~ l 7i, using any of the standard analytical approaches. 
The eigenvalue problems ffTTj) and f|T6|) can of course be studied numerically; this was done 
in Ref. 6| where the Bloch wall was found to be stable for all examined values of h. However, 
numerical solutions tend to overlook subtleties (e.g. exponentially small eigenvalues) and 
give limited insights into the structure of the configuration space. This motivates our search 
for an analytical stability proof. 

We provide such a proof in the next section. 

IV. XY MODEL IN ELLIPTIC COORDINATES 

Returning to the energy ([2]), we define elliptic coordinates on the (Reip, Im ip)-p\ane: 

xjj = B (sinh u sin v + i cosh u cos v) . (17) 

Here, u(x) > and < v(x) < 2n are continuous fields. 

The use of elliptic coordinates was pioneered by Trullinger and DeLeonardis who utilised 



9 



these in their calculation of the partition function for the MSTB model 22|. Elliptic coor- 
dinates allow the separation of variables in the effective Schrodinger equation that arises in 
their transfer- matrix approach (see also |20j). Subsequently, Ito [23J used elliptic coordinates 
to separate variables in the Hamilton- Jacob i formulation of equation (j3j) (see also 241]) . 
Transforming to the elliptic coordinates (I17p . the energy functional (j2J) acquires the form 

f\u)+g\v) 



E = 2h 



(sinh u 4~ sin v )(u x + v x ) + 



sinh u + sin v 



dx, 



where 



f(u) = B sinh u ( cosh 2 u 



A 2 



B 2 



g(v) = B sin v ( — cos 2 v 



The integrand in (Tl8l) admits a Bogomolnyi-type decomposition 

P 9/, // r Wf 9 ^ V 

E = 2h < Vj M x + — r + [V x + — r 



+ $x > dar, 



where 



and 



/i(w, f ) = sinh 2 m + sin 2 v > 



^ = -2f(u)u x - 2g(v)v x 



(18) 



(19a) 
(19b) 



(20) 



Since both terms in the square brackets in ([20]) are nonnegative, the energy is bounded from 
below: 

oo 

E > 2h$(x) . (21) 

— oo 

Evaluating the right-hand side of (121 p using (I19p . this inequality is transformed into 

B 00 
E > —{[3A 2 cosh u{x) - 5 2 cosh 3 w(x)] + [3A 2 cos v{x) - B 2 cos 3 u(ar)] } . (22) 

3 J -oo 

In terms of the elliptic coordinates, the domain walls' boundary conditions t/>(±oo) = ±iv4 
acquire the form 



//( — x } ■ air-cosh ( — 



.4 



v{— oo = 7r; 



w(+oo) = arccosh ( — ) , t>(+oo) = 0, 



(23a) 
(23b) 
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or 



u(+oo) = arccosh 




v (+00) = 2n 



(23c) 



Using these, the inequality ( |22l) becomes simply 



E[if>) > E B 



(24) 



where Eb is the energy of the Bloch wall given by j7]). 

The bound (124"1) is obviously saturated by the Bloch walls ([6]). We now show that, given 
the boundary conditions ( 1231) . the two Bloch walls are the only solutions with the minimum 
energy. The proof appeals to the Bogomolnyi equations 



which have to be satisfied by any configuration with E[ip] = E B - In order to prove the 
uniqueness, it is sufficient to demonstrate that the dynamical system (1251) has a unique 
heteroclinic trajectory connecting the point (I23al) to the point ( I23bl) and another unique 
trajectory connecting (I23al) to (I23cl) . To this end, we note that the line u = arccosh(v4/Z?) 
is an invariant manifold and that this manifold is attractive: trajectories flow towards this 
line but no trajectories can leave it. Therefore, the only trajectories connecting the points 
f l23|) have to be segments of this straight line. Letting u = arccosh ( A/ B) , equation f!25bj) 
simplifies to 



Subject to the boundary conditions v(— 00) = ir, v(oo) = and v{— 00) = ir, v(oo) = 2n, 
this equation has a unique pair of solutions 



(More precisely, these solutions are unique up to translations x — > x — x .) Inserting equa- 
tions (126!) into (TTTT) yields the right- and left-handed Bloch walls, equation (Q. Consequently, 
the Bloch walls are indeed the unique minimal energy solutions (modulo translations). This 
proves their stability — within each of the three evolution equations ([8]), Qj and (jTUl) . 




(25b) 



(25a) 



v x = —B sin v. 



sinw = ±sech(5x), cosf = tanh(Sx). 



(26) 
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V. CONCLUDING REMARKS 



We have utilised the Bogomolnyi construction to prove that the Bloch walls of equation ([2]) 
are energy-minimizing kinks. We have also shown that they are unique energy minimizers. 
Thus for all three evolution equations (JED, ©, and (ITUi) . the Bloch walls are proven to be 
stable. 

The energy (T5]) was introduced to describe the anisotropic XY model. However, the 
related evolution equations ([8]), ([9]) and ({TO]) emerge in several other areas where our results 
will also be applicable. In particular, the Ginsburg-Landau equation ([8]) appears as a generic 
amplitude equation in resonantly forced oscillatory media near the Hopf bifurcation 25) . The 
Bloch and Ising walls are often regarded as the basic building blocks for the one- and two- 
dimensional patterns arising in such media 26j. Our stability result puts this interpretation 



on a firmer ground. Next, the MSTB model (Q was studied in the context of quantum field 



theory 



271 ] . Here, the fact that the model admits a BPS bound is of fundamental importance 



as it means that it admits a natural supersymmetric extension. Finally, the parametrically 
driven NLS equation ffTUl) describes Faraday resonance in a wide shallow water channel in 



the low- viscosity limit |lOl. 28fl. Since the Ising wall has already been shown to be stable 
within the NLS equation |6(, a similar conclusion obtained now for the Bloch wall reveals 
an interesting bistability of the two solitons. This bistability should allow experimental 
realization. 



Acknowledgments 

It is a pleasure to thank Alexander Ianovsky and Dmitry Pelinovsky for useful remarks. 
S.W. was supported by a grant from the Science Faculty of the University of Cape Town. 
LB. was supported by the NRF of South Africa under grant 2053723. 



APPENDIX A: THE STABILITY-MINIMALITY CORRESPONDENCE FOR 
THE NONLINEAR SCHRODINGER DYNAMICS 

While the linearised operator coincides with the hessian of energy in the case of the 
Ginsburg-Landau and relativistic dynamics, the two operators are different in the nonlinear 
Schrodinger case. The hessian H is given by equation ffl2l . whereas the linearised operator 

12 



is J _1 7^, with J the skew-symmetric matrix (fT6|) . The multiplication by a skew-symmetric 
matrix changes the spectral properties of an operator; for example, the continuous spectrum 
of TC lies on the positive real axis whereas the continuous spectrum of J~ l Ti consists of 
pure imaginary A. Therefore it is not obvious how the energy minimality translates into the 
absence of unstable eigenvalues. 

In this appendix we provide a simple proof that the absence of negative and zero eigen- 
values of TC L is sufficient for J~ l TL not to have eigenvalues with positive real part. This 
fact is usually familiar to workers in this field; for the comprehensive treatment, including 



the necessary conditions for stability and the eigenvalue count, see 29]. 

First of all, we note that if A is an eigenvalue of the operator J~ X TL in (fiT)|) . associated 
with an eigenvector 

a(x) 



Z(.r) 

b(x) 

then —A is also an eigenvalue, associated with the eigenvector 

(a(—x) \ 
, J" (A1) 

[This follows from the fact that both for the Bloch and Ising wall, the off-diagonal elements 
of the matrix TC in ffl2|) are odd functions of x: TCwix) = H.2i(x) = — H.2i{— x).] Hence real 
eigenvalues always come in (A, —A) pairs. On the other hand, if A is a complex eigenvalue 
with an eigenvector z(x), then A* is an eigenvalue with an eigenvector z*(x). Therefore, 
complex eigenvalues appear in (A, —A, A*, —A*) quadruplets. Next, if 

Hz = \Jz, ReA^O, 

the eigenvector z satisfies an identity 

(z*, Jz) = 0, (A2) 
which follows from the self-adjointness of the operator 7i. In ( 1A2|) . we used the notation 

(zi,z 2 ) = J[ai(x)a 2 {x) +b 1 (x)b 2 {x)]dx, 
where and 6j are the (complex) components of the vector z i; i.e. 

aAx) 

z,(.r) - I | . / 1.2. 

h{x) 
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Assume now that the operator J l 7i has an eigenvalue A with Re A > 0, with the eigen- 
vector z. It is not difficult to show that the quadratic form (fl~4"l) calculated on the function 



y = Cz + C*z* + Cz + C*z*, 



(A3) 



where C and C are complex coefficients, is either sign- indefinite, or equals zero. Indeed, 
substituting flA3[) into (jl~4l) and making use of flA2[) . we get 



where z is as in (1A1I) and c.c. stands for the complex conjugate of the immediately preceding 
term. The expression ( 1A4I) is either zero or changes its sign under C — > —C. (This conclusion 
obviously remains valid if the eigenvalue A is real.) 

Thus if the form 5 2 E is positive definite on the subspace S defined by the constraint 
f fT3|) — the condition satisfied if the operator 7i does not have negative or zero eigenvalues 
other than the translation mode — the product operator J _1 7i cannot have eigenvalues with 
positive real part. 
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